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^ ! Abstract 

■^ I We present a fully relativistic formalism for describing neutrino-induced A-mediated single-pion 

t^ ' 

production from nuclei. We assess the ambiguities stemming from the A interactions. Variations 

in the cross sections of over 10% are observed, depending on whether or not magnetic-dipole 
dominance is assumed to extract the vector form factors. These uncertainties have a direct impact 
^ ' on the accuracy with which the axial-vector form factors can be extracted. Different predictions for 

a: 

C5 (Q ) induce up to 40-50% effects on the A-production cross sections. To describe the nucleus, we 
►-j^ , turn to a relativistic plane-wave impulse approximation (RPWIA) using realistic bound-state wave 

o ■ 

(y-^ I functions derived in the Hartree approximation to the a-co Walecka model. For neutrino energies 

^SJ ' larger than 1 GeV, we show that a relativistic Fermi-gas model with appropriate binding-energy 



% 



correction produces comparable results as the RPWIA which naturally includes Fermi motion, 



o 

00 . 

f~^ , nuclear-binding effects and the Pauli exclusion principle. Including A medium modifications yields 



a 20 to 25% reduction of the RPWIA cross section. The model presented in this work can be 
naturally extended to include the effect of final-state interactions in a relativistic and quantum- 
mechanical way. Guided by recent neutrino-oscillation experiments, such as MiniBooNE and K2K, 
and future efforts like MINERi^A, we present Q^, W, and various semi-inclusive distributions, both 
for a free nucleon and carbon, oxygen and iron targets. 

PACS numbers: 13.15.+g, 21.60.-n, 24.10.Jv, 25.30.Pt 



Electronic address: christophe.praet@ugent.be] 



'f Present address: Institut fiir Theoretische Physik, Universitat Giessen, Germany. 

1 



I. INTRODUCTION 



In the last few years, precision measurements of the neutrino-oscillation parameters have 
driven the interest in medium-energy neutrino physics. The MiniBooNE l| and K2K [2| 
collaborations have recently collected a wealth of neutrino data in the 1-GeV energy range 
[31], where the vast part of the strength can be attributed to quasi-elastic (QE) processes 
and A-mediated one-pion production. A thorough understanding of these cross sections is 



essential to reduce the systematic uncertainties, 
and future neutrino experiments like MINERt/A 



n turn, the high- statistics data from these 
^ and SciBooNE [sl] offer the opportunity 



to address a variety of topics related to hadronic and nuclear weak physics. 



Various theoretica 
cleon 



y,y,y,y, 



models have been developed to study one-pion production on a free nu- 



lOj. These efforts chiefly focus on studying the vector and axial- vector form 
factors that are introduced to parameterize the incomplete knowledge of the A-production 
vertex. Whereas the vector form factors can be reasonably well determined from electro- 
production data [71, Isj, the axial- vector ones remain troublesome due to the large error flags 
present in early bubble- chamber neutrino data and sizeable model dependencies in their anal- 
yses 



IG 



ll| . Besides, different theoretical calculations of the most important axial form 
factor, C^((5^), reveal highly divergent pictures 12|, [l3|, [ij] . Consequently, the Q^ evolution 
of the axial form factors and the axial one-pion mass M4 are rather poorly known. Concern- 
ing the A-decay vertex, it has been established that the traditionally-used decay couplings 
are not fully consistent with the Rarita-Schwinger field-theoretic description of the A par- 
ticle [l5|. Instead, a consistent interaction can be constructed, which couples solely to the 
physical spin-3/2 part of the A propagator [l5|. Since planned neutrino-scattering experi- 
ments aim at putting further constraints on Ma and the axial form factors, it is important 
to assess the ambiguities related to the incomplete knowledge of the A interactions. 
Modeling is made even more challenging by the fact that nuclei are employed as detectors. 
Thus, various nuclear effects need to be addressed in order to make realistic cross-section 
predictions. Traditionally, the Fermi motion of the nucleons inside the nucleus is described 
within the relativistic Fermi gas (RFC) 16|, llTJ]. Owing to its relative simplicity, the RFC 
model has been the preferred nuclear model in neutrino-event generators. Going beyond the 
?IFG, realistic bound-state wave functions can be calculated within a relativistic shell model 
171, llSj, or by adopting spectral-function approaches that extend beyond the mean- field pic- 



ture 19|, |20|, |2l|, |22| . A comparison of these models provides insight into the nuclear-model 
dependence of the computed cross sections. Another nuclear effect stems from the fact that 
the A properties are modified in a medium 23||, generally resulting in a shift of the peak 
position and a collisional broadening of the width. Finally, one must consider the final-state 
interactions (FSI) of the outgoing pion and nucleon. To study the effect of FSI, recent efforts 
have resorted to a combination of semi-classical and Monte-Carlo techniques 2J, |25|. Based 
on these results, it is clear that FSI mechanisms produce by far the largest nuclear effect on 
one-pion production computations. 

In this work, we present a fully relativistic formalism that can serve as a starting point 
to investigate A-mediated one-pion production from nuclei. Recognizing the ability of the 
new generation of experiments to measure both inclusive and semi-inclusive observables, 
we develop a framework that is geared towards a detailed study of various distributions, 
like Q^, W, energy and scattering-angle dependences. To model nuclear effects, we turn 
to the relativistic plane-wave impulse approximation, using relativistic bound-state wave 



functions that are calculated in the Hartree approximation to the a-u Wa 
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ecka model 26 1 
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29|, and 



This approach was successfully applied in QE nucleon-knockout studies 
includes the effects of Fermi motion, nuclear binding and the Pauli exclusion principle in a 
natural way. Medium modifications of the A particle are taken into account along the lines 
of Ref. f23|]. We investigate the sensitivity of the cross section to uncertainties in the A 
couplings. Then, we proceed with the nuclear-model dependence of our results, by compar- 
ing with RFG calculations. Our findings in this regard are of great importance to neutrino 
experiments that employ the RFG model in the event generators. The formalism outlined 
in this work is an ideal starting ground to implement also FSI effects in a fully relativistic 
and quantum-mechanical way. The discussion of FSI mechanisms, however, falls beyond the 
scope of the present paper. 

The paper is organized as follows. Section [III introduces the formalism for the elementary 
A-mediated one-pion production process. The third section deals with the nuclear model 
and discusses the framework for the description of neutrino-nucleus interactions. Numerical 
results are presented in section IV. In section V, we summarize our conclusions. 



II. CHARGED-CURRENT PION NEUTRINOPRODUCTION ON THE NU- 
CLEON 

A. Cross section 

For a free proton target, the charged- current (CC) process under study is 

^11 + P -^ /i" +p + 7r+. (1) 

The corresponding reactions for a free neutron are 

A+ _ n 

A+ _ _i_ ^ ' 

Isospin considerations allow one to relate the strength of the above reactions 

a{W+p ^-^ PTT+) = 9a{W+n ^ mT+) 

(3) 

where W~^ denotes the exchanged weak vector boson. In a laboratory frame of reference, 
the corresponding differential cross section is given by 30| 

1 rrii, mi d^ki tjin d^k^ d^k^^ 



d^a 



'l3E,Ei (27r)3 En (27r)3 2E^(27r)3 ^4^ 

Figure [1] defines our conventions for the kinematical variables. The target nucleon has four- 
momentum kfq^i = (mAr,0), with niN the nucleon's mass. We write ki, = {Ey.ky) for the 
incoming neutrino, ki = {Ei, ki) for the outgoing muon, k.„ = (i?^, A;^) for the outgoing pion 
and fcjv = (-Eat, k]\[) for the outgoing nucleon. The xyz coordinate system is chosen such that 
the z axis lies along the momentum transfer q, the y axis along k^x ki, and the x axis in the 
lepton-scattering plane. In Eq. (jlj), the incoming neutrino's relative velocity (3 = \ku\/ Ey is 
1. The neutrino mass rrijy will cancel with the neutrino normalization factor appearing in the 
lepton tensor. The (5-function expresses energy- momentum conservation and Xlfil^fi ^^ P 
denotes the squared invariant matrix element, appropriately averaged over initial spins and 
summed over final spins. Using the 5-function to integrate over the outgoing nucleon's three- 
momentum and the magnitude of the pion's momentum, one arrives at the fivefold cross 




FIG. 1: Kinematics for neutrino- induced charged-current one-pion production on the nucleon. 



section 



dV 



mi,mi\ki\mN\kTT\ 



dE,da,da, 2{2rfE,\EN + E4\k,r--q-K)/\K\^\ (5) 

where the soUd angles Qi and Q^^ define the direction of the outgoing muon and pion respec- 
tively. 



B. Matrix element for resonant one-pion production 



Next to the kinematic phase-space factor, Eq. ([5]) contains the squared invariant matrix 
element 

(6) 






(free) 
fi 






Here, the sum over final muon and nucleon spins is taken. Averaging over the initial nucleon's 
spin, S7v,i , leads to a factor 1/2. An explicit expression for the invariant matrix element is 
obtained by applying the Feynman rules in momentum space. Writing 



J. A free) _ GpCOS^ jp(free) ^ ^ 1" ^ 



(7) 



with Gp the Fermi constant and 6c the Cabibbo angle, one distinguishes the hadron current 



jpifree) 



■^ "^had ^~ "^V^N-i Sn)^ ^T^N^^,P^'-WNA'^y^N,ii SN,i), 



(8) 



the weak boson propagator 



S^ 



W,pa 



9paM^ 



2 ' 



Q' + M^ 



Q' 



-QpQ' 



and the lepton current 



(9) 



(10) 



Clearly, the least-known physics is contained in the vertex functions of the matrix element 
of Eq. ([8]). For the A- production vertex, we adopt the form 



rly7VA(^A,g) 



+ 



rriN "" ~^ 



m 
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m 



N 



rriN 



fyA/Q2\ ryA(Q2\ 

ig"' A - q^Y) + ^^^i^ig^'q ■ ^a - q'ki) + c^{Q')9^^ + ^^^-^g^g" 
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where a set of vector {C} 



3.. 6) and axial (C/ 



(11) 
3. .6) form factors are introduced. 



These form factors are constrained by physical principles and experimental data. Imposing 
the conserved vector current (CVC) hypothesis leads to Cq = 0. The partially-conserved 
axial current (PCAC) hypothesis, together with the pion-pole dominance assumption, yields 
the following relation between C^ and the pseudoscalar form factor Cq 



ct 



5 ,n2 I ,^2 



'AT 



Q"^ +77ll 



At Q^ = 0, the off-diagonal Goldberger-Treiman relation gives C^ 



(12) 
1.2 J9|. Furthermore, 



CVC entails that the weak vector current and the isovector part of the electromagnetic 
current are components of the same isospin current. Consequently, after extracting the elec- 



tromagnetic form factors from electroproduction data, the C, 



V 



3, 4, 5 follow immediately 



by applying the appropriate transformations in isospin space. To extract the vector form 
factors, it has been established that the magnetic-dipole (Ml) dominance of the electromag- 
netic N ^ A transition amplitude is a reasonable assumption [31|]. This Ml dominance 



leads to the conditions 



32| 



a 



V 



3 W ' 



cl 



0, 



(13) 



where W is the invariant mass, defined as 1^ = \/kj^. For C^ , a modified-dipole parame- 



terization is extracted 



Bfl 



32| 



cl 



1.95D 



V 



Qy4M^ 



(14) 



with Dv = {1 + Q^/M^)-^ the dipole function and My = 0.84 GeV. In Eq. ([HD, the 
faster-than-dipole fall-off reflects the fact that the A is a more extended object than a 
nucleon. Within this scheme, it is possible to relate all weak vector form factors to C^ . 
More recently, a direct analysis of the electroproduction helicity amplitudes from JLab and 
Mainz experiments resulted in an alternative parameterization of the weak vector form 



Is 



factors [9 1 

y _ 2.UDv y _ -1.51 y 

^3 -l + gV4M2' ^^ - 2.13 ^3' 

^ l + gV0.776M2' 

attributing a non-zero strength to the weak vector form factor C^ . The axial form factors 
are even more difficult to determine, in the sense that they are only constrained by the 
bubble-chamber neutrino data. A widely used parameterization is given by 

1.2 1 



24 






(16) 



4 ' "^3 



where M4 = 1.05 GeV. However, there still resides a great deal of uncertainty in the axial 



form factors or, equiva 



ently, in C^. The extracted axial-mass value, for example, is heavily 



model-dependent 10|, lUl]. A re-analysis [10] of ANL data within a model that includes 
background contributions next to the A-pole mechanism reveals a C^(0) value that is lower 
than the one predicted by the Goldberger-Treiman relation. This result is corroborated by 



some recent form-factor calculations, within a chiral constituent-quark (yCQ) model 12| 
and lattice QCD framework [13|, [ij] . Figure [2] compares both theoretical calculations with 
the phenomenological fit of Eq. ( TT6l) . It can be clearly seen that all three approaches exhibit 
highly divergent Q^ evolutions. 
The Rarita-Schwinger spin-3/2 propagator for the A reads 



where Ma = 1.232 GeV and F stands for the free decay width. 
A common way of describing the A decay is through the interaction Lagrangian 

C^NA = ^^^^f t(9^0)V^ + h.C, (18) 



0, 
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FIG. 2: Different results for the axial transition form factor C5 (Q^). The full line represents a 
phenomenological fit to ANL and BNL data within a A-dominance model (Eg. 1161) . The dashed 
line shows a quenched lattice result, and is parameterized as C5 (Q^) = C's (0)(1 + Q"^ / M'j<^~'^ , 
^5(0) = 0.9 and Ma = 1-5 GeV [14]. The dotted line corresponds to a xCQ result, and is taken 
from Ref. [ijl. 

where t/'^, and ip denote the spin-3/2 Rarita-Schwinger field, the pion field and the nucleon 
field respectively. The operator T is the isospin 1/2 -^ 3/2 transition operator. From (TTSll . 
one derives the simple vertex function 






m^ 



and the corresponding energy-dependent width 

V{W) - ^ '"^^ 
with 



127r m£l4/ 



(mAf + E 



iV 



l^c 






(19) 



(20) 



(21) 



Requiring that F(Ma) equals the experimentally determined value of 120 MeV, one obtains 
fnNA = 2.21. An alternative choice for the AttN interaction Lagrangian is provided by 



*-'7r 



NA 



JttNA a 



^^^"G'^a7M75Tt(9,0)^, 



(22) 



m^MA 

where G^a = Ofsipa — daipp- This form has been proposed by Pascalutsa et al. [15], who point 
out that many of the traditional couplings, like the one in Eq. (ITSl) . give rise to unwanted 
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spin-1/2 contributions to the cross section. The interaction of Eq. (122|) . however, couples 
only to the physical, spin-3/2 part of the A propagator. With the interaction Lagrangian 
of Eq. 0221) the vertex function becomes 

T'^^Akn, k^) = ^^6'^"^^A:.,„7/375fcA,-y. (23) 

Calculating the decay width from Eq. (!23l) leads to the same expression as in Eq. (l20l) . 

implying f*^^ = f^NA = 2.21. 

Combining formulas ([6]) to (ITOil . the squared invariant matrix element can be cast in the 

form 

I nAfree)i2 _ G p COS OgM^^ ^p„ . 

where the leptonic tensor is given by 

2 

with the definition 69123 = +1- Introducing the shorthand notation O'^ = ^A^^N^^,^iv^wNA^ 
one arrives at the following expression for the hadronic tensor 

H'free) = ^Tt (^{ /c^, + mN)dP{ /c^ + m^)0''^ , (26) 

where Op = 7o(C^)'^7o- 



III. CHARGED-CURRENT PION NEUTRINOPRODUCTION FROM A NU- 
CLEUS 

Turning to nuclear targets, a schematical representation of the reaction under study is 
given by 

u^ + A^fi- + {A-l) + N + 7T, (27) 

where A denotes the mass number of the target nucleus. Compared to the free-nucleon case, 
one now needs to consider the residual nucleus fc^-i = {Ea_-i, A;a-i) as an extra particle in 
the hadronic final state. Following the same line of reasoning as in section Hi Al the lab-frame 



cross section corresponding to the process of Eq. (127|) becomes 

(fa 
dEidQidET^dQ^dfiN 

nij^mi I ki I miymA-i | ^tt 1 1 ^tv | 



2{2tt)^E,\Ea-i + En + E^kN ■ (k^ - q}/\kN\ 

'ft 

A. Relativistic bound-state wave functions 



(28) 






The invariant matrix element in fl28l) carries the tag bound and involves nuclear many- 
body currents between initial and final nuclear wave functions. In medium-energy physics, 
however, one usually resorts to a number of assumptions that allow a reduction of the 
nuclear-current matrix elements to a form similar to Eq. ([S]). Here, we summarize the main 



approximations that enable this simplification and refer to Ref. [33| for the more detailed 
and analytic considerations. First, we only consider processes where the residual {A — 1) 
system is left with an excitation energy not exceeding a few tens of MeV. The major fraction 
of the transferred energy is carried by the outgoing pion and nucleon. Further, we adopt the 
impulse approximation (lA): the nuclear many-body current is replaced by a sum of one- 
body current operators, exempt from medium effects. Assuming an independent-particle 
model (IPM) for the initial and final nuclear wave functions, the hadronic current matrix 
elements can be written in the form of Eq. ([8]), whereby the initial- nucleon free Dirac spinor 



is replaced by a bound-state spinor 33j. This approach, where the outgoing nucleon and 
pion remain unaffected by the nuclear medium, is generally referred to as the relativistic 
plane- wave impulse approximation (RPWIA). 

The single-particle wave functions used in this work are determined in the Hartree approx- 
imation to the a-uj Walecka model, using the Wl parameterization for the different field 



strengths (26|. They are written as 



*.„(.! = I '"1^""-% I . (29) 

where m is the magnetic quantum number and a stands for all other quantum numbers that 
specify a single-particle orbital. In the definition of the spherical two-spinors, a generalized 
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angular momentum k is introduced. The momentum wave functions are obtained from 

^-^rn{p) = 7^^ j ^„,™(r)e-*P"W. (30) 



.(!_;) ^ /^l / 9{p)y+K,v^{p) 



with 



The result is 

9{p) = / G{r)ji(pr)dr, (32) 

Jo 

and 

/•oo _ I I _L_ \ K < Q \ 

f{p) = sgn{K) / F{r)jj{pr)dr, 1=1 ' . (33) 

>/o \^/-i, K>oy 

In ( !32|l and (!33|) . j;(x) = x ji{x) is the Ricatti-Bessel function. 

Returning to the calculation of the squared invariant matrix element in fl28|) . the following 

factor appears 

Sa{p) = ^ ^^a,m{P)^a,m{p)- (34) 

This expression, referred to as the bound-state propagator, can be cast in a form which is 
similar to the free-nucleon projection operator 3J]. One finds 



^a(p) = (^ + M„), (35) 

with the definitions 

E. = j^,jA9^P) + nP))^ (36) 

^" = (^^ {up)mp) ■ 

In other words, the hadronic tensor for scattering off a bound nucleon is readily found from 
the free-nucleon one in Eq. (126!) by making the replacement 

-^ ^ (2vr) (^ + M.). (37) 

Figure [3] shows the momentum wave functions of Eqs. (l32l) and ( l33ll for a proton belonging 
to a specified carbon shell. Owing to the small contribution of the lower wave-function 
component, the quantities Mq, and Ea are almost equal in strength. 
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FIG. 3: The left panel shows the momentum wave functions for the carbon nucleus. The full 
(dashed) line corresponds to g{p) {f{p)) for a lsi/2 proton, the dotted (dash-dotted) line represents 
dip) ifip)) foi' ^ lP3/2 proton. In the right panel, the quantities defined in Eq. (j36p are shown for 
a lp3/2-shell ^^C proton. 

B. Medium modifications of A properties 



In a nuclear environment, the A mass and width will be modified with respect to its 
free values. These medium modifications can be estimated by calculating the in-medium A 
self-energy, as was e.g. done in Ref. 23|]. The real part of the A self-energy causes a shift of 
the resonance position, whereas the imaginary part is related to the decay width. Medium 
modifications for the width result from the competition between a Pauli-blocking correc- 
tion, reducing the free decay width, and a term proportional to the imaginary part of the 
A self-energy, including various meson and baryon interaction mechanisms and, therefore, 
enhancing the free decay width. A convenient parameterization for the medium-modified 
mass and width of the A is given in Ref. [23|, in terms of the nuclear density p. For our 
purposes, we shall adopt an average nuclear density p = O.TSpo, with po the equilibrium 
density. Then, at the A peak, we calculate the following shifts 



Ma — > Ma + 30 MeV, 
F — ^ F + 40 MeV. 



(38) 
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In Ref. [35|, a similar recipe was used to accommodate medium modifications in the calcula- 
tion of ^^C(7, pn) and ^^C(7, pp) cross sections. There, the computations proved to compare 
favorably with the data in an energy regime where the reaction is dominated by A creation. 

IV. RESULTS AND DISCUSSION 

In this section, we present computations for the process 

^^l+P — ^ /i"+p + 7r+, (39) 

the strength of which can be straightforwardly related to the other channels listed in Eq. ([2]) 
by applying the isospin relations of Eq. ([3]). Throughout this work, the cross sections are 
shown per nucleon, meaning that the RPWIA results are scaled to the number of protons 
in the target nucleus. Unless otherwise stated, we use the vector form factors derived in 
the Ml-dominance model (Eqs. (fT3|) and (fT^ ). the axial form factors of Eq. (fT6l) with 
Ma = 1.05 GeV, and the traditional AttN coupling defined in Eq. flTSj) . For the RFG 
calculations, we adopt kp = 225 MeV and an average binding energy of Eb = 20 MeV. The 
latter value can be considered as a good estimate for the weighted average of the centroids 
of the single-particle strength distributions in typical even-even nuclei near the closed shells 



36|. 



A. WNA and AttN couplings 

Before discussing the nuclear effects described in section llllt we address some topics 
related to the elementary A couplings introduced in section IIIBI Figure H] appraises the 
sensitivity of the Q^ distribution to uncertainties residing in the vector and axial- vector form 
factors. The left-hand panel contrasts the widely used Ml-dominance parameterization for 
the vector form factors with a more recent fit to electroproduction helicity amplitudes, which 
extends beyond magnetic-dipole dominance p] . With the Lalakulich fit of Eq. flTSl) one finds 
cross sections which are about 10% higher than those obtained with the Ml-dominance 
form factors of Eqs. flT3|) and flT^ . The discrepancy between both parameterizations is 
significant, as vector form factors are usually regarded as well-known when they are used 
as input to extract the far less known axial form factors from neutrino-scattering data. 
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FIG. 4: Q^ evolution of the A^^-production cross sections for a lsi/2 ^^C proton and an incoming 
neutrino energy of 1 GeV. In the left panel, the full (dashed) line corresponds to the vector form- 
factor parameterization of Eqs. (J13p and (J14p (Eq. (jlSp ). The right panel studies the sensitivity of 
the cross sections to the various parameterizations for C5 (Q'^), contained in Fig. [2j 



As pointed out in section IIIBt the current situation for the axial-vector form factors is 
somewhat more dramatic. To see how uncertainties in C^(Q^) affect the cross section, we 
have performed computations with both the phenomenological result of Eq. fITB]) and the 
theoretical calculations shown in Fig. [21 The right-hand panel of Fig. H] shows what this 
implies for the Q^ distribution. Clearly, the Q^ evolution of the A-production cross section 
exhibits a strong sensitivity to the adopted C^{Q^) parameterization. Near Q^ = 0, cross 
sections using the xCQ- and QCD-model results are about 40% lower than the calculation 
with the A-dominance fit. This is almost entirely due to the difference in C^(0) values, 
which yields a ratio of (0.9)^/(1.2)^ ^ 0.56 for the dominant cross-section contribution. 
The rapid fall-off predicted by the xCQ model results in cross-section values that are much 
lower over the whole Q^ range. On the other hand, the QCD calculation foretells a less 
steep dipole dependence, leading to more strength towards higher Q^ values. The xCQ 
result for C^ halves the integrated cross section with respect to the calculation with the A- 
dominance fit. To investigate the impact of different A-decay couplings, we have computed 
ly-distributions using both the traditional coupling of Eq. (ITS]) and the Pascalutsa coupling 
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FIG. 5: Invariant-mass dependence of the A'*"+-production cross sections for a lsi/2 ^^C proton 
and an incoming neutrino energy of 1 GeV. The hadronic invariant mass is defined as T^ = 



\J{kT^ + k^Y . The full (dashed) fine uses the AvrA^ coupling of Eq. (fTHj) (Eq. 

of Eq. ( !22l) . The results are shown in Fig. [5], where it can be seen that differences between 
the two approaches are smalL Although the Pascalutsa coupling yields higher values in the 
tail of the PF-distribution, we infer an overall effect that does not exceed the 2% level. 

B. Nuclear-model effects 



In this subsection, the results of section IIVAI will be put in a more general perspective. 
To this end, we will compare neutrino-nucleus with neutrino- nucleon cross sections. Figure [6] 
shows how the total strength for the process in Eq. (1391) varies with the incoming neutrino 
energy. Under the same kinematical conditions and with similar input for the A couplings, 
our results for the elementary process compare very well with the predictions published in 
Ref . [lO| . Turning to the predictions for target nuclei. Fig. [6] shows how the elementary 
cross section is halved near threshold. For higher incoming energies, the effect dwindles to 
20% at E^ = 800 MeV and 8% at E^ = 2 GeV. Most strikingly, the RFG calculations are 
in good to excellent agreement with both the carbon and iron RPWIA results. The only 
discernable feature of Fig. E] is that the iron curve exceeds the carbon and RFG ones by 
roughly 15% just beyond threshold. This can be understood after recognizing that the iron 
result is largely due to outer-shell protons, which are less bound than the corresponding 
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FIG. 6: Total cross sections per nucleon for i^fj, + p -^ n +p + 7r+. The full line represents the 

elementary process, for scattering from a free proton. The dash-dotted line stands for the RFG 

calculations, whereas the dashed (dotted) line corresponds to scattering from a carbon (iron) target 

nucleus. The right panel focusses on the threshold region. 

carbon ones. Clearly, the nuclear-target cross sections are very sensitive to binding-energy 
differences at lower incoming energies. These effects, however, vanish at higher neutrino 
energies and fall to a 1% correction level at E,y = 1 GeV. As a matter of fact, at sufficiently 
high energies RFG calculations with a well-chosen binding-energy correction are almost 
indiscernible from the corresponding RPWIA results. These findings are more detailedly 
assessed in Figs. U\ M and M Figures [7] and [S] compare RFG and RPWIA computations. 
The former considers scattering from a carbon target at E^, = 800 MeV, which corresponds 
to the mean energy of the neutrino beam used by the MiniBooNE experiment. As can be 
appreciated from Fig. [3, the RFG and RPWIA models produce almost identical results. In 
Fig. [HI we present the ratio of RFG to carbon RPWIA results for the twofold cross section 
(P a / dT^^d cos 6* , where T^ is the outgoing pion's kinetic energy and 6* its scattering angle 
relative to the neutrino-beam direction. Apart from the threshold region, where numerical 
instabilities induce large fluctuations, it is observed that differences between the RFG and 
RPWIA result do not exceed the 5% level over the whole (T^r, 6*) range. In addition, the 
largest deviations occur where the cross section has hardly any strength. Consequently, 
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v^+p -^ fj- +p + 7r+. A carbon target and an incoming neutrino energy of 1 GeV are considered. 

upon integrating over T^ and ^*, we find that the total RFG cross section exceeds the 
RPWIA one by about 2%. Figure M compares the cross section for a carbon nucleus with 
the one for an iron nucleus at E„ = 1.5 GeV. Although the total strength, integrated over 
the outgoing muon energy Ei, is the same for both nuclei, it is interesting to note that the 
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FIG. 9: The i^^ + p ^ fi + p + Tr"*" cross section per nucleon as a function of the lepton energy 

El for Eu = 1.5 GeV. The fuh hne represents the elementary process, whereas the dashed (dotted) 

line refers to scattering from carbon (iron). 

iron distribution is shifted with respect to the carbon cross section. Again, this reflects 
the fact that a carbon proton requires, on average, more energy than an iron proton to 
be knocked out of the nucleus, leaving therefore less energy for the outgoing muon. To 
assess the influence of medium modifications to the mass and width of the A, we have 
computed the ly-distribution according to the prescription given in Eq. ( l38l) . In Fig. [10, 
we contrast the elementary cross section with RPWIA calculations for carbon. Compared 
to the free case, the RPWIA result that does not include A medium modifications is seen 
to be heavily suppressed at the A pole. The nuclear binding brings about a broadening of 
the py-distribution, mainly reallocating strength to lower invariant masses. When medium 
modifications are taken into account, the A pole is shifted towards higher W values, by an 
amount that corresponds to the mass shift in Eq. (1381) . Compared to the plain RPWIA 
case, the peak is again suppressed and broadened, owing to the increased medium-modified 
width. On the whole, we observe a 25% reduction of the RPWIA cross section when A 
medium modifications are included. 
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incoming neutrino energy of 1 GeV. The full (dashed) hue represents the calculation for a free 
(bound) proton. The dotted curve adds the effect of A medium modifications to the RPWIA 
result. 

C. Results under MiniBooNE and K2K kinematics 

In view of recent results presented by the MiniBooNE and K2K collaborations [3|, we 
conclude this section with some computations for the specific neutrino energies and target 
nuclei employed by these experiments. From an experimental viewpoint, the most accessible 
distributions are the ones with respect to outgoing-muon variables. Fig. [H] depicts an 
RPWIA calculation for a two-fold differential cross section against the outgoing-muon energy 
and scattering angle with respect to the neutrino beam. The incoming neutrino energy is 
fixed at 800 MeV, corresponding to MiniBooNE's mean beam energy. Since MiniBooNE has 
carbon as target material, this calculation was performed on a carbon nucleus. The result 
shown in Fig. [H] can be integrated over 6i or Ei to yield the one-fold cross sections displayed 
in Fig. [121 Relative to the free cross section, the angular distribution for a carbon target is 
evenly reduced by about 20%. In general, the outgoing muon prefers a forward direction, 
although a minor shift seems to take place between the free and the bound case. This effect 
relates to the change in the muon-energy distribution, depicted in the right-hand panel of 
Fig. [T2j Indeed, for scattering off bound protons, one observes a shift of the Ei distribution 
towards lower values. Recognizing the correlation between high muon energies and forward 
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FIG. 11: Cross section per nucleon for f^+p -^ /x +p + 7r+ against outgoing-muon energy and 
scattering angle. The incoming neutrino energy is 800 MeV, the target nucleus is carbon. 

scattering angles, as can be appreciated in Fig. [TTl the bound case will correspondingly 
yield a larger number of events at slightly higher scattering angles. We also note that the 
RPWIA result fades out sooner than the elementary cross section, because a certain amount 
of energy is needed to knock the carbon proton out of its shell. Planned experiments like 
MINERz/A endeavor to have a good energy resolution for both the muon and the hadronic 
final state. The ability to detect the outgoing pion or nucleon or even both would allow 
a detailed study of different nuclear effects. In Figs. [13] and [H] we present cross sections 
versus the pion kinetic energy T^^ and pion scattering angle relative to the beam direction 
6*. This time, we adopted K2K settings, namely an oxygen target hit by neutrinos with an 
energy of 1.3 GeV. From the left-hand panel of Fig. [HI one infers that, within the RPWIA 
model, the outgoing pion prefereably leaves the nucleus along the beam direction. As for 
the kinetic-energy distribution, we observe a comparable reduction and shift of the strength 
as in the muon-energy distribution. 

V. CONCLUSION AND OUTLOOK 



We have developed a relativistic framework to study A-mediated one-pion production 
from nuclei at medium energies. The proposed formalism offers great flexibility in calcu- 
lating various observables both for the free process and for scattering from nuclear targets. 
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FIG. 12: Cross sections per nucleon for i/^ +p -^ fi + p + vr^ , for 800 MeV neutrinos scattering 

from a carbon target. The left (right) panel shows the cross section as a function of the outgoing- 

muon scattering angle (energy). Each of the panels contrasts the elementary cross section (full 

line) with the RPWIA result, without A medium modifications (dashed line). 

Motivated by operational and planned experiments, we have conducted a systematic study 
by addressing the impact of A-coupling ambiguities on the Q^ and W distributions. Cross 
sections are found to vary by as much as 10% depending on whether or not the Ml-dominance 
assumption is used to extract the vector form factors. This is very significant, as the ex- 
tracted value for the axial mass Ma depends heavily on the model applied in the analysis of 
the neutrino-scattering data and, therefore, on a reliable input for the vector form factors. 
Uncertainties in the dominant axial form factor, C^((5^), have a dramatic effect on the 
A-production cross sections. At low Q^, a 25% reduction of the off-diagonal Goldberger- 
Treiman value C^(0) = 1.2 leads to cross sections that are smaller by 40%. In the case of a 
Q^ dependence that is steeper than a modified- dip ole form, the effect increases to almost 50% 
over the whole Q^ range. In the W distribution, we observe 2%-level deviations between the 
traditional A-decay coupling choice and a consistent one, which effects a decoupling from the 
spin- 1/2 terms. To investigate the influence of nuclear effects, we have computed RPWIA 
neutrino-nucleus cross sections for carbon, oxygen and iron nuclei. We have briefly touched 
on the topic of A medium modifications. Using a prescription that gives good results in 
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FIG. 13: Cross section per nucleon for f^+p — > fi +p + 7r+ against outgoing-pion kinetic energy 
and scattering angle. The incoming neutrino energy is 1.3 GeV, the target nucleus is oxygen. 



photo-induced two-nucleon knockout and electron-scattering studies, we infer a 20-25% sup- 
pression of the RPWIA cross sections due to medium effects. The nuclear responses are 
very sensitive to binding-energy differences at lower neutrino energies. From E^ = 1 GeV 
onwards, the cross sections per nucleon for different nuclear targets are seen to agree at the 
1% level. To assess the nuclear-model uncertainty in our description of A-mediated one-pion 
production, we have also contrasted the RPWIA results with calculations performed within 
an RFG model with a well-considered binding-energy correction. At 1-GeV neutrino ener- 
gies, differences between one- and two-fold distributions computed within both models do 
not exceed the 5% level. The agreement is better for total cross sections, where deviations 
between the RFG and RPWIA model dwindle to 1-2%. Hence, for sufficiently high incom- 
ing neutrino energies, the influence of Fermi motion, nuclear binding and the Paul exclusion 
principle can be well described by adopting an RFG model with binding-energy correction. 
The RFG model, however, just as the RPWIA approach, falls short in implementing FSI 
and nuclear correlations of the short and long-range type. Contrary to the RFG, the model 
proposed in this work has the important advantage that it can serve as a starting point for 
a relativistic and quantum-mechanical study of FSI mechanisms. As a matter of fact, the 
inclusion of FSI for the ejected pions and nucleons is currently under study. To this end, 
we closely follow the lines of Ref. [33|, where use is made of a relativistic Glauber model for 
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FIG. 14: Cross sections per nucleon for y^+p — > fi" +p + 7r^, for 1.3 GeV neutrinos scattering 
from an oxygen target. The left (right) panel shows the cross section as a function of the outgoing- 
pion scattering angle (kinetic energy). Each of the panels contrasts the elementary cross section 
(full line) with the RPWIA result, without A medium modifications (dashed line). 

fast ejectiles and an optical-potential approach for lower ejectile energies. 
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